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Low-frequency K^ = 0"*" states in deformed neutron-rich nuclei are investigated by means of 
the quasiparticle-random-phase approximation based on the Hartree-Fock-Bogohubov formalism in 
the coordinate space. We have obtained the very strongly collective K'^ — 0"*" modes not only 
in neutron-rich Mg isotope but also in Cr and Fe isotopes in A'^ = 40 region, where the onset of 
nuclear deformation has been discussed. It is found that the spatially extended structure of neutron 
quasiparticle wave functions around the Fermi level brings about a striking enhancement of the 
transition strengths. It is also found that the fluctuation of the pairing field plays an important role 
in generating coherence among two-quasiparticle excitations of neutron. 

PACS numbers: 21.10.Re; 21.60.Ev: 21.60.Jz 



I. INTRODUCTION 

Physics of nuclei far from /3-stability has been one of 
the main current subjects of nuclear physics. Quest for 
new kinds of collective motion in exotic unstable nuclei 
is one of the most interesting issues in nuclear structure 
physics and has been actively studied both experimen- 
tally and theoretically [1]. Because low-lying collective 
excitations are sensitive to the shell structure around the 
Fermi level, one can expect unique excitation modes to 
emerge associated with the new spatial structures such as 
neutron skins and novel shell structures which generate 
new regions of deformation. 

Breaking of the spherical magic number iV = 20, 28 
and striking enhancement of B{E2; — > 2^) in Mg iso- 
topes towards drip line are under lively discussions in 
connection with onset of the quadrupole deformation and 
continuum coupling 0, i, i, i, i, 0, S, i, M EH, IH] ■ Fur- 
thermore in heavier mass region, collectivities and break- 
ing of the iV = 40 sub-shell have been recently discussed 
for neutron-rich Cr and Fe isotopes [H, [M [l^, [H, d ■ 
Although many investigations based on the random 
phase approximation (RPA) or the Quasiparticle RPA 
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(QRPA) including the pairing correlations [20, i21l 

IMM y M l2dM l29l Im, [12, S l3^, m 

on multipole responses of neutron-rich nuclei have been 
done, they are largely restricted to spherical systems. 

Recently, low-lying RPA modes in deformed neutron- 
rich nuclei have been investi gat ed by several groups [H, 
SO, SH, El S H EE m, \M\M, \M ■ These calculations, 
however, do not take into account the pairing correlation, 
or rely on the BCS pairing (except for Ref. [i^ and our 
very recent investigation [49|), which is inappropriate for 
describing the pairing correlation in drip line nuclei due 
to the unphysical nucleon gas problem [50|. 

In Ref. [51], we constructed a new computer code 
that carries out the deformed QRPA calculation based 



on the coordinate-space Hartree-Fock-Bogoliubov (HFB) 
formalism, and investigated quadrupole excitations in 
neutron-rich Mg isotopes close to the drip line around 
N = 28. We obtained low-lying K'^ = 0+ and 2+ 
quadrupole excitation modes in 36,38,40]y[g found 
that the K"^ = O"*" mode in *°Mg was generated by co- 
herent superposition of neutron two-quasiparticle (2qp) 
excitations in loosely bound and resonant states. Among 
them, the neutron 2qp excitations of up-sloping oblate 
levels and down-sloping prolate levels play the major 
role. A deformed gap is formed at = 28 around 
/32 = 0.3 [1^] due to the crossings between the up-sloping 
i/[303]7/2 level and the down-sloping i/[310]l/2 level, and 
this deformed closed shell approximately corresponds to 
the (/7/2)~^(P3/2)^ configuration in the spherical shell 
model representation (see Fig. [1]). It was also found 
that the dynamical pairing correlation, i.e., fluctuation 
of the pairing field played a crucial role in generating 
coherence among neutron 2qp excitations. 

The single-particle levels which are (2j -I- 1)— fold de- 
generate in the spherical potential are split by the nuclear 
deformation to 2-fold degenerate levels with the single- 
particle energy to first order of (32 as [l^ 



3»^~j(j + l) 



(1) 



In a prolately deformed potential (/32 > 0), a level with 
small value of f2 (z— component of the angular momen- 
tum j) has a negative slope. We call this single-particle 
level a down-sloping prolate level in this paper, and refer 
a positive slope level with large 11 as a up-sloping oblate 
level. 

Because the deformed shell gap is, in many cases, made 
by the up-sloping and the down-sloping levels , we in- 
vestigate in this paper detailed and systematic properties 
of low-frequency K'^ ~ 0+ modes in deformed neutron- 
rich Mg, Cr and Fe isotopes, and discuss generic features 
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of the soft K'^ = modes uniquely appeared in de- 
formed neutron-rich nuclei. 

The paper is organized as follows. In the next section, 
we present our model of the deformed HFB and QRPA. 
Numerical results and discussion are given in § IIIII Fi- 
nally, we summarize the paper in § IIVI 



II. MODEL 

We briefly summarize our approach (see Ref. [Slj for 
details). In order to discuss simultaneously effects of nu- 
clear deformation and pairing correlations including the 
continuum, we solve the HFB equations [50, 52, 53|| 



-ih-{ra 



VT,a('^cr)^ 



'a 



(2) 



directly in the cylindrical-coordinates assuming axial and 
reflection symmetries. Here t — v (neutron) and tt (pro- 
ton), and r = {p,z,(f>). For the mean-field Hamiltonian 
h, we employ the deformed Woods-Saxon potential with 
the parameters used in Ref. for '^^Mg and those in 
Ref. [39J for neutron-rich Cr and Fe isotopes. The pair- 
ing field is treated self-consistently by using the density- 
dependent contact interaction 54|, ISq. 



Vpp{r,r') = Vo 



I -Pa 



g^^(r) 
Qo 



5{r-r'). (3) 



with Qq = 0.16 fm~'^. Here Q^^{r) denotes the isoscalar 
density and Pa the spin exchange operator. Because the 
time-reversal symmetry and the reflection symmetry with 
respect to the x — y plane are assumed, we have only to 
solve for positive 51 and positive z. We use the lattice 
mesh size Ap = Az — 0.8 fm and the box boundary 
condition at /Omax = 10.0 fm and Zmax = 12.8 fm. The 
quasiparticle energy is cut off at 50 MeV and the quasi- 
particle states up to O'^ = 15/2=*= are included. 

Using the quasiparticle basis obtained by solving the 
HFB equation ([2]), we solve the QRPA equation in the 
matrix formulation ^5^ 



7^ 



9-,s 



1 

-1 



/q/3 



. . (4) 

The residual interaction in the particle-particle (p-p) 
channel appearing in the QRPA matrices A and B is the 
density-dependent contact interaction ([3]). On the other 
hand, for the residual interaction in the particle-hole (p- 
h) channel, we use the Skyrme-type interaction 



Vphir,r') = 



io(l + XoPa) + j{l + X3Pa)g'^{r) 




FIG. 1: Single-particle energies in the deformed WS po- 
tential for neutrons in ''"Mg, plotted as functions of the 
quadrupole deformation parameter /32. Solid and dotted 
lines denote positive- and negative-parity levels, respectively. 
Single-particle levels are labeled with the asymptotic quan- 
tum numbers [NnsA]fl. 



with to = -1100 MeV-fm^, = 16000 MeY-im^,xo = 
0.5, and 2:3 = 1.0 [53|. Because the deformed Wood- 
Saxon potential is used for the mean field, we renor- 
malize the residual interaction in the p-h channel by 
multiplying a factor fph to get the spurious = 1+ 
mode (representing the rotational mode) at zero energy 
{vph fph ■ Vph)- We cut the 2qp space at Ea + Ep < 30 
MeV due to the excessively demanding computer mem- 
ory for the model space consistent with that adopted in 
the HFB calculation. Accordingly, we need another fac- 
tor fpp for the p-p channel. We determine this factor 
such that the spurious = 0+ mode associated with 
the particle number fiuctuation appears at zero energy 
(vpp — > fpp ■ Vpp). The dimension of the QRPA equa- 
tion (HI) is 2512 for ^'^Mg and 3804 for '^^Pe. Within 
this model space, the energy-weighted sum rule for the 
quadrupole K"^ — 0+ excitation is satisfied 92.9% for 
34Mg and 97.0% for ^^Fe. 



III. RESULTS AND DISCUSSION 



Let us now discuss the low-frequency = 0+ modes 
in ■^''Mg in ij llll Al and in neutron-rich Cr and Fe isotopes 



with TV ~ 40 in § IIITB] 



<5(r-r'), 
(5) 



Figure [T] shows the single-particle energy diagram for 
the WS potential as a function of deformation parameter 
f32- At around f32 — 0.5, the deformed shell gap is formed 
at = 22 due to the crossings between the up-sloping 
i/[202]3/2 and the down-sloping z^[330]l/2 and i/[321]3/2 
levels. In Fig. ^ we show the low-lying excitation spec- 
trum below 3 MeV in '^'^Mg. Here excitation energies are 



3 



2 



34 

Mg 

(4+) 


/ 4+ 


■>+ 

1 

2+ 

2+ 

K''=l* 

0+ 


2+ 


— r 




0+ 


0* 

TV Belyaev 




Exp. 


Cal. 





FIG. 2: Excitation energy spectrum in Mg obtained by the 



QRPA calculation using P2 ~ 0.4 and Vo 
and available experimental data 4^. 
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FIG. 3: QRPA strength distributions for the A"" = 0+ 
isoscalar quadrupole p-h excitations (left) and the monopole- 
and quadrupole-pair excitations (middle and right) in '^''Mg, 
calculated using fph = 0.74 and fpp — 1.37. For comparison, 
unperturbed 2qp transition strengths are shown in the lower 
panels. The arrow indicates the neutron threshold energy 
Eth ~ 4.55 MeV for Iqp continuum. 



O2 states; 



E{I, K) = hujBPA + 



2Jtv 



+ (6) 



with the vibrational frequencies, wrpa, and the Thouless- 
Valatin moment of inertia, J7tV; calculated by the QRPA. 
The pairing strength Vq for ^''Mg is determined in or- 
der to reproduce the experimental excitation energy of 
i?(2]^) = 660 keV In this figure, we also show the 
excitation energies evaluated by using the Belyaev mo- 
ment of inertia, jTsciyaov, where the residual interaction 
is turned off. We found that JTtv becomes about 11% 
larger than J/Boiyacv due to the time-odd component in 
the residual interactions of ([3]) and (O- The K"^ = 2+ 
mode at 2.01 MeV is dominantly generated by the proton 
p-h excitations as that in 36,38,40]y[g ^j^ggg g^pg ppQ_ 
ton dominated 7— vibrational modes. The negative par- 
ity K'^ = 0~ state at 2.61 MeV is mainly generated by 
the lowest neutron 2qp excitation of i^[202]3/2(X)j/[321]3/2 
with the excitation energy of 2.74 MeV. Its collectivity 
is thus rather weak. 

The left panel of Fig. [3] shows the K"^ — 0+ isoscalar 
quadrupole transition strengths, and the unperturbed 
2qp transition strengths are shown in the lower panel. 
The intrinsic transition strength to the lowest K'^ = 0+ 
excited state is about 43 Weisskopf units (1 W.u. ~ 6.5 
fm'' for ■^'^Mg). This state is constructed by the coher- 
ent superposition of 2qp excitations of the up-sloping 
(j/[202]3/2)2 and the down-sloping (i/[330]l/2)2 and the 
(i/[321]3/2)^ levels, together with the high-lying 2huj ex- 
citation of protons. 

Generation mechanism of the soft = 0+ mode in 
deformed neutron-rich nuclei is understood essentially by 
the schematic two- level model in Ref. [s^ . They consider 
the case where only two AA components are present in the 
wave functions both of the ground 0+ and of the excited 



10+) = ^L=|AiAi) + ^=L=|A2A2) (7a) 



^ |AiAi) + ^=l=|A2A2). (7b) 



= 



+ 62 ' 



The transition matrix element for the quadrupole opera- 
tor is then 



(oilQ2o|o+) = 



2a6 



2 ' a2 -f 62 



[(Ai|Q2o|Ai)-(A2|Q2o|A2)] (8) 



and it is proportional to the difference in the quadrupole 
moments of the individual orbitals composing the 0"*" 
states. In the case that the quadrupole moments of the 
orbitals have opposite signs to each other, this matrix 
element becomes large. This situation is realized in the 
level crossing region, where the up- and down-sloping or- 
bitals exist. As the number of components increases in 
the QRPA calculations, the wave function becomes more 
complicated, and the transition matrix element reads (see 
Appendix B) 



(i^^ = o+|g2o|o) 



a/3 af3 



(9) 

In the QRPA employing the separable interaction, for 
the major component of the QRPA eigenmode associated 
with the pairing fluctuation, the phase of the QRPA am- 
plitudes, faf} and gaf}, are opposite between particle- like 
levels and hole-like levels (due to the {u^ — v'^) factor; 
see Appendix J of Ref. [60i|). In the level crossing re- 
gion, the particle-like level is the up-sloping oblate level, 
and the hole-like level is the down-sloping prolate level. 
Because the transition matrix elements for the QRPA 
eigenmodes are determined by the sum of products of 
the QRPA amplitudes and the individual 2qp transition 
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TABLE I: QRPA amplitudes for the = 0+ state at 1.37 MeV in ^''Mg. This mode has the proton strength B{E2) = 
15.9e^fm^, the neutron strength B{Q^2) — 162fm'', and the isoscalar strength B{Q^^2) = 279fm^, and the sum of backward- 
going amplitudes ^ |3a/3p = 0.141. The single-quasiparticle levels are labeled with the asymptotic quantum numbers [A^nsAJfi. 
Only components with /^^ — g^p > 0.01 are listed. 
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the lowest /C" = 0+ state at 1.37 MeV. Sohd and dotted lines 
denote the isoscalar and neutron transitions, respectively. 



matrix elements in Eq. ([5]), they contribute coherently 
for the quadrupole transition. 

In the present case, the main components generating 
the low-lying K'' = 0+ state are [202] 3/2 (g) [202]3/2, 
[321]3/2 (g [321]3/2 and [330]l/2 (g [330] 1/2, which are 
the up-sloping and down-sloping levels with opposite 
quadrupole moments and the QRPA amplitudes have op- 
posite signs as shown in Table HI Therefore, the strength 
for the quadrupole transition becomes enhanced. These 
three 2qp excitations contribute about 50 % of the to- 
tal transition strength. Furthermore, many other 2qp 
excitations coherently participate to generate the lowest 
0+ state, which brings about further enhancement of the 
transition strength. 

Figure [3] shows the partial sum of the transition matrix 
element for the K'' = 0+ mode J^e^+e^^e^^^ 
the isoscalar and neutron excitations. We see that many 
2qp excitations construct this state with coherence; each 
2qp transition matrix element {M^^^) has a same sign. 
If the QRPA mode were a destructive state or a single 
2qp excitation, the sum of the transition matrix element 
becomes zero in total or shows step function. 

From this figure, we see another interesting feature: 
Soft K'^ — mode in "^^Mg is generated not only 
by neutron 2qp excitations around the Fermi level with 
the excitation energy around 5 MeV (three 2qp ex- 
citations discussed above) but also proton excitations 



around 20 MeV. These proton excitations are 2huj exci- 
tations of giant resonance; 7r[211]3/2 7r[431]3/2 (19.5 
MeV), 7r[101]3/2 7r[321]3/2 (17.3 MeV), 7r[110]l/2 
7r[330]l/2 (19.6 MeV) and 7r[220]l/2 ^ 7r[440]l/2 (20.1 
MeV). Due to this coupling with giant resonance of pro- 
ton, the proton transition strength B{E2) is larger than 
that in ^OMg; B{E2) = 15.9e2fm4 in ^'^Mg and SAe^fm'^ 
in ^°Mg [5l|. Hindrance of the proton excitation in drip- 
line nuclei is due to the extreme spatial extension of neu- 
tron wave functions and decoupling with proton wave 
functions [sH ]. 

The discussion for enhancement of the quadrupole 
transition strength to the low-lying K'^ = 0"*" mode is 
also applicable to the pair transition strength. In the 
middle and right panels of Fig. [31 we show two kinds of 
strength for neutron pair transitions; one is the monopole 
pair transition and another the quadrupole pair transi- 
tion, 

Pjo= J dr^^irr)^j^ir i), (10a) 
pL = J drr^Y2o{f )^Hr i)- (10b) 

For the quadrupole pair transitions, we see a promi- 
nent peak at 1.37 MeV, whereas no peak is seen at 
this energy for the monopole pair transitions. The ori- 
gin of this contrasting behavior is understood as fol- 
lows: The quadrupole (monopole) pairing matrix ele- 
ments for individual 2qp excitations are in opposite phase 
(in phase) between the down-sloping prolate levels and 
the up-sloping oblate levels. As the transition matrix el- 
ements for the QRPA eigenmodes are determined by the 
similar way for the p-h transition matrix elements (See 
Eq. (B4)), they contribute coherently and destructively 
for the quadrupole and monopole pair transitions, respec- 
tively. As a consequence, the quadrupole-pair transition 
strength to the lowest I-C^ = Q'^ collective excited state is 
greatly enhanced, whereas the monopole-pair transition 
strength to this state almost vanishes. 

Therefore, in a deformed system where the up- and 
down-sloping orbitals exist near the Fermi level and the 
pairing fluctuation becomes important, we can expect the 
emergence of the low-lying mode which has extremely 
enhanced strengths both for the quadrupole p-h and for 
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FIG. 5: Transition densities to the — 0^ state at 1.37 
MeV. The particle-hole transition densities of protons and 
neutrons (upper panels), the particle-particle and hole-hole 
transition densities of neutrons (lower panels). Solid and 
dashed lines indicate positive and negative transition densi- 
ties, and the contour lines are plotted at intervals of 3 x 10"* 
fm~'^. The thick dotted lines denote the proton and neutron 
half densities, 0.041 and 0.049 fm~^, respectively. 



the quadrupole p-p (pair) transition. 

The low-lying K'^ — 0+ modes possibly emerge both 
in stable and in unstable deformed nuclei, and the unique 
feature in neutron-rich nuclei is that they are mainly gen- 
erated by neutrons, whose wave functions have spatially 
extended structure due to the shallow Fermi level. In or- 
der to see the unique spatial structure, we show in Fig. [5] 
three kinds of transition densities; 

Sg^^Hp, z) = (A| i;l{p, z, a)MP, ^, ^M, (Ha) 

(J 

PPP('')(P, z) = {Xli^lip, z, \)^l{p, z, i)|0), (lib) 
P"^('')(p,^) = (A|V'.(p,^,i)V'.(p,2,T)|0), (He) 

where they are called the p-h transition density, the p-p 
(particle-pair) transition density, and the h-h (hole-pair) 
transition density, respectively [23. The upper panels 
show a typical shape vibration of protons and neutrons 
along the symmetry axis (z— axis). We can clearly see the 
spatial extension of the neutron wave functions; the exci- 
tation takes place in the skin region around 7-8 fm. This 
spatial extension of the neutron wave functions brings 
about the enhancement of the transition strength. An- 
other unique feature of the low-lying K'^ — O"*" mode is 
that the particle-pair and hole-pair transition densities 
have an appreciable amplitude and are almost compara- 
ble to the p-h transition density. 
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FIG. 6: Same as Fig.[T]but in '^"'Cr. 



B. Cr and Fe isotopes around A'^ = 40 

The deformation region near the shell closure gives a 
favorable situation for emergence of the soft = O'^ 
modes discussed in the previous subsection because many 
orbitals with opposite quadrupole moments are crossing 
around the Fermi level. Since the low-lying K'" = 0"^ 
mode is quite sensitive to the pairing correlation and de- 
formation, we investigate in the present subsection the 
low-frequency = 0^ excitations in neutron-rich Cr 
and Fe isotopes in = 40 region, and show generic fea- 
tures of the soft K'^ = 0+ modes in deformed neutron- 
rich nuclei. 

In a prolately deformation region around (32 — 0.3, 
deformed shell gaps associated with different /32 values 
appear at various neutron numbers due to the cross- 
ings between the down-sloping z/[440]l/2, i/[431]3/2 and 
i^[431]l/2 levels and the up-sloping z/[301]3/2, i/[301]l/2 
and iy[303]5/2 levels (see Fig.^. 

According to the Skyrme HFB calculations [6l|, [63 |. 
the ground state of the neutron-rich Cr and Fe isotopes 
with A^ ~ 40 has a finite deformation around P2 — 0.3. 

The top panels in Fig. [7] show the isoscalar quadrupole 
transition strengths for the K"^ = 0+ excitations. We 
can see a prominent peak at around 1 MeV in all nuclei 
under consideration. All low-lying states obtained here 
have extremely enhanced transition strengths (note that 
1 W. u. is 14.6-17.1 fm'i for 62-70cj. ^nd 15.2-17.8 fm"* 
for ^^~^^Fe). These excitation modes are generated by 
coherent superposition of neutron 2qp excitations among 
the up- and the down-sloping levels as in '^^Mg. 

Furthermore, we can see an interesting feature from 
this systematic calculation; the transition strengths be- 
come enhanced at A^ = 38 and 46, and are symmet- 
ric at N = 42. This feature can be understood by 
the role of the Q'^ = 1/2+ levels. In ^^Cr and ^''Fe, 
2qp excitations of the up-sloping (i^ [301] 3/2)^ and the 
down-sloping (i/[440]l/2)^ have largest contributions to 
the lowest K'^ — 0+ state. The transition strength of 
the (i^[440]l/2)2 excitation is 86 fm" (5.9 W.u.) in ^^Ci 
and 64 fm'* (4.2 W.u.) in ^"^Fe, whereas that of the 
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FIG. 7: Isoscalar quadrupole (top), the monopole-pair (middle), and the quadrupole-pair (bottom) transition strengths for the 
K'' = 0+ excitations in '^^-^Cr and ^'^''^'^Fe obtained by the QRPA calculations using P2 = 0.3 and Vo = -420 MeV-fm^. The 
self-consistent factors are /ph=0.72-0.74 and fpp = 1.43 for Cr isotopes, and /ph=0.74-0.79 and /pp=1.39-1.42 for Fe isotopes. 
The arrows denote the neutron threshold energies. 



(j^[301]3/2)2 is 9.6 fm'^ and 8.6 fm'*, respectively. This 
indicates the [440] 1/2 level has a spatially extended 
structure (in ^^Cr, the root-mean-square radius of the 
i'[440]l/2 level is 5.1 fm, whereas the total neutron r.m.s. 
radius is 4.2 fm). As the neutron number increases, the 
transition strength decreases as a consequence of the de- 
creasing contribution of the (i/[440]l/2)^ excitation. 

As the neutron number approaches N — 46, the 
[431] 1/2 level becomes close to the Fermi surface (see 
Fig. In ^"Cr and ^^Fe, 2qp excitations of the up- 
sloping (i^[303]5/2)2 and the down-sloping (i^[431]l/2)2 



have largest contributions to the lowest = 0+ state. 
The transition strength of the (z^[431]l/2)^ excitation is 
115 fm" (6.7 W.u.) in ™Cr and 95 fm'' (5.4 W.u.) in ^^Fe. 
This 2qp transition strength becomes large because of the 
shallower Fermi level and the spatially extended structure 
of the quasiparticle wave function of the z/ [431] 1/2 level. 

At = 42, both of the [440] 1/2 and the z/ [431] 1/2 
levels are located far from the Fermi level, and 2qp ex- 
citations of the (z/[301]l/2)2 and the (z^[422]5/2)2 have 
main contributions to the lowest = 0+ state in ^^Cr 
and ^^Fe. These quadrupole transition matrix elements 
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have opposite signs, but these transition strengths are 
less than 1 W.u. 

In Fig. [3 we show also the strengths for the monopole- 
and quadrupole-pair transitions. As in '^'^Mg, these low- 
lying = 0^ modes have enhanced strengths for 
the quadrupole pair transition, whereas disappearing 
strengths for the monopole pair transition. This system- 
atic calculation shows an importance of the dynamical 
pairing: Fluctuation of the pairing field, which is de- 
formed as well as the mean field, is crucial for generating 
the collective K'^ = 0"*" modes. 



IV. SUMMARY 

We have studied low-frequency K"^ — 0+ modes in 
neutron-rich nuclei, taking account of the effects of nu- 
clear deformation, pairing correlation and continuum 
coupling simultaneously. New type of this excitation 
mode is generated by coherent superposition of neutron 
2qp excitations near the Fermi level whose wave func- 
tions have spatially extended structure. It is found that 
the dynamical pairing correlation, i.e., the pairing vibra- 
tion enhances its collectivity. 

In the spherical neutron-rich nuclei, the effect of the 
dynamical pairing has been investigated in detail in 
Refs. [1^ [23, ll^l- We have shown in this paper the 
importance of the dynamical pairing in neutron-rich de- 
formed systems. In a deformed system where the up- and 
down-sloping orbitals exist near the Fermi level, one ob- 
tains the low-lying mode possessing extremely enhanced 
strengths both for the quadrupole p-h transition and for 
the quadrupole p-p (pair) transition induced by the pair- 
ing fluctuation. 

We found that the coupling between the pairing vibra- 
tion and the neutron-skin vibration brings forth the soft 
K"^ = 0+ mode in deformed Mg region with N = 22. 
Furthermore, it was shown that emergence of low- lying 
K"^ — O"*" modes is not restricted in the neutron-rich Mg 
isotopes. As an example, in neutron-rich Cr and Fe iso- 
topes around iV = 40, we showed that the coherent cou- 
pling between the pairing vibration and the /5— vibration 
of the neutron skin brings about the striking enhance- 
ment of the strengths for the quadrupole p-h and the 
quadrupole pair transition. 
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APPENDIX A: MATRIX ELEMENTS FOR 
ONE-BODY OPERATORS 

Let us first consider matrix elements for particle-hole 
type one-body operators 

(a6|6(,"^)|HFB), (Al) 

where 

0(^") =11/ drdr' 5„,„,5{r-r')0'';^^\r)i;\r' a')ij{ra), 

(A2) 

and the HFB ground state and the 2qp excited states 

/3,|HFB)=0, (A3a) 
\ah) = plpl\R¥B) (A3b) 

are described by the quasiparticle operators. These op- 
erators are defined by the generalized Bogoliubov trans- 
formation 

i^\ra) ^ ^i^k{ra)pl + (plk{r<j)$k, (A4a) 
fe 

Mra) = V^lkir^)^k + '/'2,fc(ra)/3^, (A4b) 
fe 

where 

tpkira) ^ -2(T(pkir -cr). (A5) 

In Appendices, we omit the subscript r for simplicity. 
The 2qp transition matrix elements are calculated as 

(a6|d^"'^|HFB) 
= (HFB|/3fa/3,6^^^|HFB) 

kk' era' 

= J drO^;^''\r){-^,Ar l)^2Ar^)+ViAr^W2Ar i) 
+ Vi,b{r i)<P2,a(r T) - Vifi{r T)¥'2,a(r- i)} 

(A6) 

using the quasiparticle wave functions, and we employed 
the Wick's theorem. In the cylindrical coordinate repre- 
sentation, we can rewrite as 

(a6|0^"^|HFB> = 
27r(5K,o„+Oi, JpdpdzO^^''\p, z) 

X WiAp^ ^) ])'^2,b{p, Z, i) - <y3l,a(p, Z, i)<^2,6(p, Z, ]) 

- 'PiAP^ ^> T)'^2.q(P, Z, i) + (filAP, Z, i)<P2,a(p, Z, t)}, 

(A7) 
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where 



(A8) 



Next we consider the pair creation operators consisting 
of nucleons with opposite direction of spins 

6^") = [ drdr'5{r-r')0^;^''\r)i;\r' T)^^(r- i). (A9) 



The matrix elements read 

(a6|6^"^|HFB) 
= (HFB|/3b/3aO(^"^|HFB) 

= ^ /" drdr'5{r-r')0'^P\r) 
fcfc' 



'/:'i.fc(r-'f)(/?i.fc'(rI)(HFB|/3,/3,/3T/3T,|HFB) 
drO(,"")(r){-^i,„(r i)^i,b(r t) + fi^r i)^i,a(r T)} 



=27r(5A',f2„+nb y pdpdzO''^'^^ {p, z) 



X {(/?l,a(p, 2, DtPl-fjlP, i) - fl,a{p, Z, l)ipi,b{p, Z, t)}, 

(AlO) 



where 



OT\p,z)^OT\r)e 



(All) 



APPENDIX B: QRPA TRANSITION MATRIX 
ELEMENTS 



In terms of the nucleon annihilation and creation op- 
erators in the coordinate representation, the quadrupole 
operator is represented as 

Q2K = Y.j drr-'Y2,-K{f)i^\ra)iXrcT). (Bl) 



The intrinsic matrix elements {\\Q2k\^) of the 

quadrupole operator between the excited state |A) and 
the ground state |0) are given by 

(Alg^KlO) = Q2KUfa\ + g'ab) = E ^4kL (b2) 

ab ab 

where 

Q2K,ab = (a6|Q2i^|HFB) (B3) 



calculated by using the quasiparticle wave- functions as 
(|X7)). 

We can furthermore calculate the intrinsic matrix ele- 
ment of the pair creation operator as 



ab 



ab 



(B4a) 



ah 



(B4b) 



where the monopole- and quadrupole-pair creation oper- 
ators are defined by Eq. (fTO|) . The matrix element of 2qp 
excitation Pg^o^j^, P^l^J^^,, P^'^l^ and Pax'^if, are given by 



(lAlOp and by the similar expression for O 



(vv) 
K ■ 



APPENDIX C: CALCULATION OF THE 
MOMENT OF INERTIA 

We estimate moments of inertia using the spurious so- 
lution of the RPA equation [l3| 



A B 

B* A* 



J X 



0, 



a/375 \ ^/ 7(5 
A B\ f e \ ^flj_ (J: 



(Cla) 
, (Clb) 



where is the angular momentum operator consisting of 
the orbital part Ix and the spin part s^^ and its conjugate 
operator 8. 

The Thouless-Valatin moment of inertia JTtv is deter- 
mined thorough the orthonormal condition 



(C2) 



as 



Jtv = 2te E {J^)IM - B)-l^s{J^U. (C3) 

Turning off the residual interaction, we obtain the expres- 
sion for the Inglis-Belyaev moment of inertia jTBoiyaev 



a/3 



I -^a^ la/3 



(C4) 
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